Hierarchy Bloch Equations for the 

Reduced Statistical Density 
Operators in Canonical and Grand 
Canonical Ensembles 

Liqiang Wei 

Institute for Theoretical Atomic, Molecular and Optical Physics 
Harvard University, Cambridge, MA 02318 

Chiachung Sun 
Institute of Theoretical Chemistry, Jilin University 
Changchun, Jilin 130023 P. R. China 

February 2, 2008 



Abstract 

Starting from Bloch equation for a canonical ensemble, we deduce 
a set of hierarchy equations for the reduced statistical density operator 
for an identical many-body system with two-body interaction. They 
provide a law according to which the reduced density operator varies 
in temperature. By definition of the reduced density operator in Fock 
space for a grand canonical ensemble, we also obtain the analogous 
Bloch equation and the corresponding hierarchy reduced equations 
for the identical interacting many-body system. We discuss their pos- 
sible solutions and applications. 
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1 Introduction 



The fundamental problem in quantum statistical mechanics is to determine 
the density matrix of a statistical ensemble. Once the density matrix is 
known, the thermodynamic properties can be calculated from the correspond- 
ing microscopic mechanical observables [1-4]. The basic equation of motion 
for the statistical density matrix is von Neumann-Liouville equation. For an 
equilibrium statistical ensemble, its solution takes the following special form, 

P = Y^Pi\^i >< A\, (1) 

i 

where pi is the thermal probability distribution, and \ipi > is the energy 
eigenstate of the system. The summation over % is for all the eigenstates 
consistent with symmetry constraints. For a canonical ensemble, for exam- 
ple, its relative thermal probability distribution and density operator are as 
follows, 

P t = eM-PEi), (2) 

and 

D N = Y,eM-PEiM >< A\ = exp(-(3H N ), (3) 

i 

where j3 is the reciprocal of the product of Boltzmann's constant kb and the 
absolute temperature T, and H N is the Hamiltonian of an iV-body system. 
The trace of Eq. (3) gives the partition function of the canonical ensemble, 

Tr(D N ) = ^exp(-/3^) = Z(J3,V,N). (4) 

i 

From their definition in Eqs. (1) or (4), we observe that, to determine 
the density matrix or the partition function, there are generally two steps 
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involved. In the first step, we have to solve an energy eigenequation for the 
interacting many-body system to get its energy eigenstates and eigenvalues. 
In the second one, we need to sum over all the microscopic states. Obviously, 
these are both a very difficult or prohibitive task for the system considered, 
especially when the quantum effects are important and can not be neglected. 
Nevertheless, some other approaches or approximations have been developed, 
including the one called the equation of motion method, to address these 
issues [5,6]. Furthermore, for an identical interacting many-body system, it 
is known that its full density matrix conveys much more information than is 
necessary, and its reduced density matrix suffices to describe the states of a 
whole system. In analogous to that for the full density matrix, we attempt 
here to develop similar equations for the reduced density matrix so that they 
provide a kind of laws that the reduced density matrix obeys, or from which, 
the reduced density matrix can be directly determined. 

This paper is organized as follows. In the next section, we derive the 
reduced Bloch equations for the reduced statistical density operators for a 
canonical ensemble. They are the equations of motion according to which 
the reduced density matrix varies in temperature. In Section 3, we define 
the reduced density matrix for a grand canonical ensemble, and obtain the 
analogous Bloch equation. We also extend the ideas developed in Section 2 
on the reduced Bloch equation from the canonical ensemble to the case of a 
grand canonical ensemble. In Section 4, we summarize our results and discuss 
the possible routes for solving the equations as well as their applications. 
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2 Bloch Equation for the Reduced Statistical 
Density Operators in Canonical Ensemble 

Consider a canonical ensemble for an interacting iV-body system whose den- 
sity matrix is shown in Eq. (3). Differentiation of both sides of Eqs. (2) or 
(3) with respect to (5 leads immediately to the differential equation, 

- ^D N = H N D N , (5) 

which is called the Bloch equation [5,6]. It provides a law according to which 
the canonical density operator varies in temperature. It is equivalent to the 
definition (3). Nevertheless, it offers an alternative route to determine the 
density matrix. If the Hamiltonian of the system is known, we can directly 
solve the exact or approximate density matrix without recourse to the two 
procedures as we outline in the Introduction Section. A noticeable feature of 
Eq. (5) is that after its substitution of f3 by it/h, it yields an equation that 
is identical in form with the time-dependent Schrodinger equation. There- 
fore, some developed techniques for solving the time- dependent Schrodinger 
equation can be borrowed to obtain the solutions for the Bloch equation (5). 

In practical applications, moreover, the interactions among the identical 
particles such as electrons in atoms, molecules, and solids are often two- 
body forces. The iVth-order density matrix contains much more information 
than is required. In order to calculate the thermodynamic properties of the 
system, we need at most the second-order reduced density matrix and not 
the full iVth-order one. That is to say, the pth-order reduced density matrix 
can describe the N(> p > 2) particle states. For this purpose, we try to find 
a method for directly solving the pth-order reduced density matrix instead 
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of the one for a full iVth-order density matrix. 

The pth-order reduced density matrix of an iV-particle system is generally 
defined by [7-11] 

D> = L* If (D»), (6) 

where L P N is the contraction operator [12,13]. Its trace gives the partition 
function defined before 



Tr(D p ) = Z(P,V,N). 



(7) 



Assume that the Hamiltonian of the system considered is made up of two 
contributions, 

N N 

(8) 



N N 

H N = J2Hi)+Y,9(hj), 

i=l i<j 



where h(i) and g(i, j) are one- and two-particle operators, respectively. Rewrite 
the Hn into 



where 



N p N N 

h n = h[+ J2 M*) + E E g(hj)+ E rtJ), 

j=p+l i=l j=p+l i<j(i>p+l) 



(9) 



(10) 



i=i 



Then substituting the Hamiltonian Eq. (9) into the Bloch equation (5) and 
applying on both sides the operator L P N , we get the following equations for 
the reduced density operators, 



= HlD p +{N-p)Ll +l [h{p+l)D p+1 ] + {N-p)Ll +l 
+ ( N : P )L p p+2 \g(p + l,p + 2)D^ 



J2g(i,p+l)D p+1 



i=l 



+ 



(11) 
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We call them the reduced Bloch equations for the canonical ensemble. They 
relates the D p , D p+1 , and D p+2 in a certain way, and can be seen as a set 
of hierarchy. From the general definition for the reduced density matrix (6), 
they can also be regarded as an independent equation for D p+2 . Similar pro- 
cedures have lead to the hierarchy equations for the reduced density matrix of 
the pure states before from the Schrodinger equation in a density matrix form 
[14-17]. It follows from the idea in deriving the Bogoliubov-Born-Green-Yvon 
(BBGKY) hierarchy in classical statistical mechanics [4]. 

3 Bloch Equation for the Reduced Statisti- 
cal Density Operators in Grand Canonical 
Ensemble 

To provide more avenues for studying the behaviors of statistical density 
matrices, the above description can also be extended to the case of a grand 
canonical ensemble. In iV-particle Hilbert space V N , the relative grand 
canonical distribution reads 



P i (N)=exp[-P{E i - l jiN)], 



(12) 



and the density matrix is 



D G (N) 



exp[-P(H - »N)], 



exp(-pB), 



(13) 



where 



H = H-nN, 



(14) 
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is called the grand Hamiltonian on V^. Again, upon differentiating both 
sides of Eq. (13) with respect to (3, we have 

-±D G (N)=HD G (N), (15) 

which has the same form as Eq. (5). We call it the Bloch equation for the 
grand canonical ensemble on iV-particle Hilbert space V N . 

For the grand canonical ensemble, the number of particles N fluctuates. 
Its density matrix should be defined in the entire Fock space: F = J2n ®V N . 
That is, it should be written as the direct sum of every density matrix D G (N) 
belonging to the iV-particle space V N , 

oo 

D G = £ <BD G (N), (16) 

7V=0 

and its trace in Fock space is the grand partition function, 

oo 

Tr(D G ) = ]T Y,<xp[-P(Ei-iiN)] = E(l3,ii,V). (17) 

N=0 i 

For the same consideration, we define the corresponding pth-order reduced 
density matrix as follows, 

D P G = t ® ( ^ ) L P N [D G (N)], (18) 

N=p \ V J 

( N \ 

where the binomial coefficient I I is added to reflect the symmetric prop- 
erty of the density matrix in terms of the permutation of particles. Obviously, 

D G = Z((3,t,,V), (19) 

and 

Tr(D G ) = /( N )\Dl (20) 



Assume Hamiltonian H in H takes the same form as in Eq. (8). Rewrite H 
into 

N P_ N_ _N_ 

(21) 



N p N N 

H — Hf + £ fc(i) + E E + E 9(i,j), 

i=P+l i=l j=P+l ?<j(i>p+l) 



where 



and 



h(i) = h(i) — (j,, 



m = jrh(i) + irg(ij). 

i=l i<j 



(22) 



(23) 



Substituting the Eq. (21) into Eq. (15) and applying on both sides the 
operator J2n= p © ^ ^ ^ ^at> we finally obtain the equations for the reduced 



density matrix 

~^D p G = mD P G + (p+lK +1 [Hp + l)D p G +1 ] + (p+l)Ll +1 



v_ 

i=l 



52g(i,P + l)I% 



p+i 



+ 



+ 



P + 2 



T p 



g(p+l,p + 2)D p G 



p+2 



(24) 



which are the reduced Bloch equation on Fock space for the grand canonical 
ensemble. They also represents a law accoding to which the grand reduced 
density matrix changes in temperature. Similar notes can also be made as in 
the case for the canonical ensemble. 



4 Summary and Discussions 

In this paper, within the framework of equilibrium quantum statistical me- 
chanics, we have derived the hierarchy equations of the density operators 
for both canonical and grand canonical ensembles. They are the equations 
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of motion in terms of the variation of temperature that the reduced density 
operator satisfies. They provide the route for a direct determination of the 
reduced density operator for the statistical ensembles. 

In the first place, the Eqs. (11) or (24) can be solved in an approximate 
way. Since they are a set of hierarchy equations, we have to relate the D' p+1 
and the D p in a certain way to decouple them. The most obvious choice is to 
take an orbital approximation. For an identical fermion system, this means 
to express a pth-order reduced density matrix as a Grassmann product of p 
first-order reduced density matrices as follows 

D p+i = DPAD^D 

= D 1 A D 1 A ... A D 1 J (D°) P . (25) 
p+i 

Similar scheme has been applied for solving the reduced Schrodinger equa- 
tion for the reduced density matrix to obtain the usual Hartree-Fock equation 
[14-17]. We expect that it will lead to more general equation beyond that for 
zero temperature when it is applied to the Eq. (24) [21]. Consequently, this 
will provide a possibility for studying the microscopic structure of high tem- 
perature species or macroscopic molecules at different temperatures [19-22]. 
In addition, it will serve as a staring point to develop more general schemes 
for studying such important issues as the interplay between the microscopic 
structure and macroscopic thermodynamic properties [18]. This also includes 
the pressure effects on the electronic structure for solids [23]. Furthermore, 
Eqs. (11) or (24) can also be solved exactly for some model systems. One 
example is to consider a pairing Hamiltonian, and solve the equations in 
an occupation number representation [25-27]. The pairing Hamiltonian is 
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a useful model, for example, in the study of the nuclear structure and the 
microscopic states of superconductors. We hope that we can obtain its exact 
reduced density matrices with different orders and therefore investigate the 
corresponding thermodynamic properties [28]. Work along these two lines 
can be carried out. 
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